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Abstract 

The current research on credit risk is primarily focused on modeling default probabilities. 
Recovery rates are often treated as an afterthought; they are modeled independently, in 
many cases they are even assumed constant. This is despite of their pronounced effect on 
the tail of the loss distribution. Here, we take a step back, historically, and start again from 
the Merton model, where defaults and recoveries are both determined by an underlying 
process. Hence, they are intrinsically connected. For the diffusion process, we can derive 
the functional relation between expected recovery rate and default probability. This relation 
depends on a single parameter only. In Monte Carlo simulations we find that the same 
functional dependence also holds for jump-diffusion and GARCH processes. We discuss 
how to incorporate this structural recovery rate into reduced-form models, in order to restore 
essential structural information which is usually neglected in the reduced-form approach. 

Keywords: Credit risk. Loss distribution. Value at Risk, Expected Tail Loss, Stochastic 
processes 
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Introduction 



An accurate description of portfolio credit risk is of vital interest for any financial insti- 
tution. It is also a prerequisite for realistic ratings of structured credit derivative products. 
Furthermore, it is a crucial aspect in banking regulations. We can distinguish two conceptu- 
ally different approaches to credit risk modelhng: structural and reduced-form approaches. 
The structural models go back to Black and Scholes (1973) and Merton (1974). The Merton 
model assumes that a company has a certain amount of zero-coupon debt which becomes 
due at a fixed maturity date. The market value of the company is modelled by a stochastic 
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process. A possible default and the associated recovery rate are determined directly from this 
market value at maturity. In the reduced-form approach default probabilities and recovery 
rates are described independently by stochastic models. The aim is to describe the depen- 
dence of these quantities on common (macroeconomic) covariates or risk factors. For some 



well known reduced-form model approach 


es see, e.g., Jarrow and TurnbuU 


(1995 


), 


Jarrow 


et al. 


(1997 


), 


Duffie and Singleton 


(1999 


), 


Hull and White 


(2000 


) and 


Schonbucher 


(2003). 



First Passage Models constitute a third approach — which is usually regarded as structural, 
but is better described as a mixed or pseudo-structural approach. They were first intro- 



duced by Black and Cox (1976). As in the Merton model, the market value of a company 
is modelled as a stochastic process. Default occurs as soon as the market value falls below 
a certain threshold. In contrast to the Merton model, default can occur at any time. In 
this approach, the recovery rate is not determined by the underlying process for the market 
value. Instead recovery rates are modelled independently, for example, by a reduced-form 
approach (see e.g., Asvanunt and Staal ( 2009a|[b )). In some cases recovery rates are even 
assumed constant, for instance, in Giesecke (2004). The independent modelling of default 



and recovery rates leads to a serious underestimation of large losses. 

In this paper we take a step back, historically, and revisit the Merton model, where 
defaults and recoveries are both determined by an underlying process. Hence, they are in- 
trinsically connected. For a correlated diffusion process the Merton model has been treated 
analytically, e.g. 



m 



Bluhm et al. 


(2002 


) and in 


Giesecke 


(2004 



lation we can also derive the functional relation between expected recovery rate and default 
probability. This relation depends on a single parameter only. In Monte Carlo simulations we 
find that the same functional dependence also holds for other processes like jump-diffusion 
and GARCH. We discuss how to incorporate this relation into reduced-form models, in order 
to restore essential structural information which is usually neglected in the reduced-form 
approach. 

The paper is organized as follows. We give a short introduction to the Merton model 
in Section [T| In Section [2] we treat the diffusion case analytically and compare the results 
to Monte-Carlo simulations. In Sections |3] and |4] we extend the Merton model to jump- 
diffusion and GARCH processes. We discuss the applicability of the structural recovery rate 
beyond the Merton model in Section |6| We summarize our findings in Section |6| 

1. Merton model 

The Merton model assumes that a company k has a certain amount of zero-coupon 
debt; this debt has the face value Fk and will become due at maturity time T. The company 
defaults if the value of its assets at time T is less than the face value, i.e., if Vk(T) < Fk. 
The recovery rate then reads Rk = Vk (T) / Fk and the loss given default is 



Vk{T) 
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We denote the loss given default with an asterisk to distinguish it from the loss including 
non-default events. The individual loss can be expressed as 



...(.-lf))e(.-^), 

where B is the Heaviside function. In the Merton model, defaults and losses — and hence 
also recoveries — are directly determined by the asset value at maturity. Therefore, the 
stochastic modelling of the market value Vk{t) of a company allows to assess its credit risk. 
Let pvi^{Vk{T)) be the probability density function (pdf) of the market value at maturity. 
Then the default probability is given by 

Fk 

PD,k = JpVk{Vk{T))dVk{T) (3) 



and the expected recovery rate can be calculated as 

Fk 

(Rk) = ^ PVkiVkiT)) m{T) . (4) 

-T D,fc ^ -Tfc 

Let us now consider a portfolio of K credit contracts, where the market value of each com- 
pany k is correlated to one or more covariates. Conditioned on the values of the covariates we 
obtain different values for ^ and (Rk)- In fact, we find a functional dependence between 
default probability and recovery rate. This is in stark contrast to modelling approaches 
which assume an independence of these quantities. In section |2] we derive this functional 
dependence analytically for the diffusion case. We compare the result with Monte Carlo 
simulations. In section [3] we extend the model by including jump terms in the random pro- 
cess. In section |4] we consider a GARCH process. In both cases we find the same functional 
dependence between default and recovery rates. 



2. Correlated diffusion 

A. Analytical discussion 

In the diffusion case we can easily derive all results analytically. To keep the notation 
simple we consider a homogeneous portfolio of size K with the same parameters for each 
asset process, and with the same face value, Fk = F, and initial market value, Vfc(O) = Vq. 
We model the time evolution of the market value of a single company A; by a stochastic 
differential equation of the form 

^ = fidt + ^/^adW^ + VT^adWk . (5) 

This is a correlated diffusion process with a deterministic term ndt and a linearly correlated 
diffusion. The parameters of this process are the drift constant /i, the volatility a and the 
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correlation coefficient c. The Wiener processes dWk and dW^ describe the idiosyncratic and 
the market fluctuations, respectively. 

For discrete time increments At = T/N, where the time to maturity T is divided into 
steps, we arrive at the discrete formulation of Eq. (|5]). The market value of company k at 
maturity can be written as 

N 

VkiT) = n (l + /"^^ + V^crri^,tVAt + - caek,tVAt) . (6) 
t=i 

We define the market return as the average return of all single companies k over the 
time horizon up to maturity. 
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k=l 
K N 



7^ E n (l + /^^^ + (^^m,tVAt + Vl-c aek,tVAt) - 1 . (8) 
^ k=it=i 

For A' — )► oo we can express the average over k as the expectation value for e^ f. Due to the 
independence of t for different k and t, we can write 

N 

X^ + l = l[(l+ fiAt + v^ar/^,tv^ + Vl^(j{ek,t)VAt) , (9) 
t=i 

with 

{ek,t) = . (10) 

Thus, expression (|9| simplifies to 

N 

t=i 



exp f 51 (-'- + crVra,t^/At 
\t=i 



^ exp\^\^fi- —jT + aVcAtJ2Vm,tj ■ (11) 

In the last step of the calculation we expanded the logarithm up to ffist order in At. The 
random variables ?7m,i are standard normal distributed. Therefore the variable 

/ ca2\„ 1 ^ 



In (X^ + 1) = - _ j T + VcT^ E Vm,t (12) 



is normal distributed with mean fiT — \ca'^T and variance ca^T. This implies a shifted 
log-normal distribution for the market return itself. 



For a single company k we can write 



In 



t=i 



^In (l + fxAt + ^/cari^^ty/^ + Vl - caek,tVAt) 



■ I 1 ^ 

-r + (7V(i-c)r^ V, 



ln(X^ + l) 



:i - c)a 



(14) 
(15) 



Conditioned on a fixed value for the market return Xm, all variables Vfc(T) are independent 
and lnVfc(T)/Vo is normal distributed with mean In (Xm + 1) — |(1 — c)o"^T and variance 
(1 — c)a^T. Since we consider a homogeneous portfolio, we will omit the index k in the 
following. This allows for a clearer notation. The probability density function for the 
market value V{T) is given by the following log-normal distribution 



Pv{V{T)) 



V{T)^2n{l-c)a^T 



exp 



/ 



V 



(in m _ In (x^ + 1) + 1 (1 _ c)a^Ty 
2(1 - c)a^T 



(16) 



We obtain the individual default probability by integrating this pdf from up to the face 
value F, 



r 

Pb(X^) = jpv{V{T))dV{T) 



' l n4-ln(X^ + l) + |(l-c)a^r ^ 



(17) 



where $ is the cumulative standard normal distrubution. The expectation value for the 
individual loss given default L* — 1 — V{T)/F can be calculated as 



'ln^-ln(X^ + l) + |(l-c)(72r^ 



(18) 



- exp In {X^ + 1) - In 



ln^-ln(X^ + l)-|(l-c)(72r' 



;i-c)a2T 



The expected recovery rate is simply 

{R{X^)) = 1 - {L*{X^)) . 



(19) 



And the portfolio loss for a homogeneous portfolio, or simply the average loss, is obtained 
as 

{L{X^)) = P^{X^){L*{X^)) . (20) 
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For the sake of clarity, we introduce the function 

A(X^) = ln^-ln(X^ + l) (21) 

and the compound parameter 

B = ^(1 - c)a^T . (22) 
Now, the expressions for Pi^^X^) and {R{X^)) simphfy to 

and 

The relation between default probability and expected recovery rate depends only on B and 
is parametrized by A{X^). Thus, the parameter B can be calibrated to default and recovery 
rate data. Furthermore, by inverting expression (23), we can express A in terms of Pd, 

A = B^-\Pn)-^BK (25) 
This leads us to the functional dependence of recovery rate and default probability, 

(i?(PD)) = -^exp(-i?$-i(PD) + ^5')$($-'(PD)-fi) . (26) 

It is worth repeating that this functional relation depends on a single parameter B only. 
In Figure [l] we demonstrate this parameter dependence of Equation (26). We plot (i?(PD)) 
for a set of parameter values. Larger values for B lead to an overall decrease of recovery 
rates. In addition the dependence on Pd becomes steeper. From Equation (26) we obtain 
the functional relation of portfolio losses and default probabilities, 

(L(Pd)) = Pd - exp (-B<^-\Pj,) + 1-bA $ (^-^(Pd) - b) . (27) 



The parameter dependence of this relation is also shown in Figure [TJ 

In order to calculate Value at Risk and Expected Tail Loss we make use of the substitution 

\pL{L)dL\ = \pxjX^)dXj . (28) 

Here we work with a simplified notation where L refers to {L{X^)). With the cumulative 
distribution function 

X 

FxJx) = J pxjX^)dX^ (29) 

— oo 
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Figure 1: Parameter dependence of the functional relations {R{Pb)) (left plot) and (L(Pd)) (right plot). 
The parameter values are B = 0.1, 0.2, 0.3, 0.4 and 0.5 (top to bottom in left plot, and bottom to top in 
right plot). 



we can express the Value at Risk VaRo, as 

VaR, = L(F^i(l-«)) , (30) 



where the function L refers to the Xm-dependence given in Equation (20) and the value 
Fxl^(l — a) is the (1 — a)-quantile of the market return. The Expected Tail Loss is calculated 
as 

1 

ETL, = - / LpL{L)dL = - f L{X^)pxjX^)dX^. (31) 
a J a J 

VaRa -oo 

We will now compare these analytic results with Monte Carlo simulations for a finite 
portfolio size. 

B. Monte Carlo results 

In the Monte Carlo simulation we consider the stochastic process in Equation ^ for 
discrete time increments. The simulation is run with an inner loop and an outer loop. In 
the inner loop we simulate K = 500 different realizations of ek,t for a single realization of 
the market fluctuations r/m,* with t = 1,...,N. The inner loop can be interpreted as a 
homogeneous portfolio of size K, or simply as an average over the idiosyncratic part of the 
process. In each run of the inner loop, we calculate the market return X^, the number of 
defaults N-o^X^^) and the expected recovery rate {R{X^)). The market return Xm is defined 
as the average return at maturity, 



^m = ^ER7^-l • (32) 




For sufficiently large K the idiosyncratic terms average out and the market return X^ is 
solely defined by the realization of rj^^f. This is the reason why we use the market return as 
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a parameter for the other observables. The number of defauhs A^D(-^m) simply counts how 
many times the condition Vfc(T) < F is fulfilled. We can estimate the default probability as 

Pd(X^) ~ N^{Xr^)/K . (33) 

We obtain the portfolio loss as the average of individual losses in Equation ([2]), 

{L{X^)) = ^j2Lk. (34) 

k=l 

Using the relation 

(L(X^)) = Pd(X^) (1 - {R{X^))) (35) 
we can estimate the expected recovery rate as 

Here, we assume that the number of defaults is strictly non-zero, which is justified for large 
portfolio size K. The outer loop runs over 10^ realizations of the market terms, where we 
obtain different values for the market return and, consequently, the number of defaults 
A^D(-^m), the default probability PD(-^m) and the expected recovery rate {R{X^)). 

In this paper we do not discuss the parameter dependence of the models or aim at 
calibrating them to a given portfolio. Instead we only present the results for a single set of 
parameters with economically sensible values. As correlation coefficient we choose c = 0.5. 
The parameters for the diffusion process are n = 0.05 and a = 0.15. The initial market 
value is set to Vq = 100, the face value of the zero-coupon bonds is F = 75 with maturity 
time T = 1. 

Figure |2] shows the dependence of recovery rates {R{X^)) on default probabilities Pj:,{X^) 
in one plot, and the dependence of portfolio losses (L(Xjn)) on default probabilities PD(-'^m) 
in a second plot. In both cases we observe a very good agreement between the Monte Carlo 



simulations and the analytical results in Equations (26) and (27), respectively. Deviations 
from the average values are smaller for higher default probabilities. Low default probabilities 
correspond to positive or small negative market returns X^, see Figure |3} In this case 
defaults and recoveries are mostly influenced by the idiosyncratic part of the process and 
the correlation to the market plays only a minor role. This is why we observe a much broader 
range of recovery rates for low default probabilities. 

When only limited historical data on defaults and recoveries is available, it might appear 
reasonable to assume a constant recovery rate. This corresponds to a linear dependence 
of the portfolio loss on default probability. For small default probabilities this is a good 
approximation of the results shown in Figure [2j However, for larger default probabilities 
a constant recovery rate model severly underestimates portfolio losses. The combination 
of high default probabilities and low recovery rates strongy influences the tail of the loss 
distribution. 
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Figure 2: Dependence on default probabilities Pq for the diffusion process. The left plot shows the 
dependence of recovery rates {R{X^)) on default probabilities PuiX^). The right plot shows the dependence 
of portfolio losses {L{X^)) on default probabilities PD(Xm)- The Monte Carlo results for the diffusion 
process (grey symbols) are compared to the analytical results (solid line). Black star symbols indicate local 
average values of the simulation results. The parameter B — -^7(1 — c)a^T is determined from the simulation 
parameters. 



In Figure |3] we demonstrate that the dependence of defauh probabihties and portfoho 



losses on the market return is well described by Equations (17) and (20), respectively. Given 
the functional relation between portfolio loss L and market return we can transform the 
pdf of market returns into the pdf of portfolio losses, 

Pl{L)= |^/(^^)| P^n.(^m)- (37) 

A comparison to the Monte Carlo result for the loss distribution is given in Figure |4j We 
observe an excellent agreement, even for extremely large portfolio losses. Finally, we present 
the results for Expected Loss, Value at Risk and Expected Tail Loss in Table |Bj The latter 
two are calculated for the 0.99-quantile. The analytical values describe the simulation results 
very well. 

Since the diffusion process was used both in the Monte Carlo simulations and in the 
analytical derivation, the good agreement of the respective results comes as no surprise. In 
the following sections we will explore two different processes in the Monte Carlo simulations 
and compare these with the analjd;ical results for the diffusion. 



3. Correlated jump— diffusion 

A. Analytical discussion 

We extend the diffusion model of the previous section by adding two jump terms to 
the stochastic process, dJk for idiosyncratic jumps and dJm for jumps which affect the 
entire market. The jump terms are not contained in Merton's original model; they ensure 
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Figure 3: Dependence on market returns for the diffusion process. The left plot shows the dependence 

of default probabilities PD(Xm) on market returns The right plot shows the dependence of portfolio 

losses {L{X^)) on market returns X^. The Monte Carlo results for the diffusion process (grey symbols) are 
compared to the analytical results (solid line). The parameter B = — c)a'^T is determined from the 
simulation parameters. 
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Figure 4: Probability density function Pl{L) for the portfolio losses for the diffusion process. The black 
solid line corresponds to the simulation results. The grey dashed line shows the analjd:ical results, where 
the pdf of market returns is transformed. 





siiiiulatioii 


aiiah'tic-al result 


EL 


7.38 ■ 10-^ 


7.35 - 10-^ 


VaRo.99 


1.30-10-2 


1.29 - 10-2 


ETLo.99 


2.38-10-2 


2.37 - 10-2 



Table 1: Expected Loss, Value at Risk and Expected Tail Loss for the diffusion process. 
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that the default probabihty does not vanish as the time to maturity becomes very short. 



For apphcations of jump-diffusion processes in credit risk modelhng see, e.g., Zhou (2001) 



Schafer et aL (2007) and Kiesel and Scherer (2007). The stochastic differential equation for 



the market value reads 
dl^ 



^dt + a/c(T dPFm + dJm + Vl - cadWk + dJk 



{31 



This is a correlated jump-diffusion process. We model the jumps by a Poisson process with 
intensity A. We recall that in such a process the probability function for the event to occur 
n times between zero and time t is given by 



P 



Poisson 



it) 



{xty 



exp (—At) 



(39) 



The size A of the jump, measured in units of the current market value Vfc(t), is a random 
variable with a distribution which we have to specify. Jumps can be positive or negative. 
The largest possible negative jump is 100% of the current market value. Based on this 
information, a possible distribution of the jump size A is a shifted lognormal distribution, 
A + 1 ~ LN(/ij, (Tj), with mean /xj and standard deviation aj. 

Without the jump term, the distribution of the market value Vk{t) is log-normal. The 
jumps render the tails of the distribution fatter. The parameters of the jump process can be 
adjusted in order to match the tail behavior of a given empirical time series of the market 
value. Here, we use the same parameters for idiosyncratic and market wide jumps. 



Let us now consider the stochastic process in Eq. (38) for discrete time increments At 



T/N, where the time to maturity T is divided into N steps. Then the market value of 
company k at maturity is given by 



N 

Vk{T) = ^0 n (l + y^^t + Vc(Jr]^,t^t + dJ^^t + ^l-caek,t^t + dJk^t 



t=i 



(40) 



Akin to the diffusion case, we find for the market return of the jump-diffusion 



1 

K 



K N 



= ^ ^ n (l + /"^^ + V^(^Vm,tVAt + dJ„,,t + VT^aek,tVAt + dJ^ J - 1 . (41) 

k=l t=l 



For K ^ oo we can express the average over k as expectation values for ek,t and dJ^,*- Since 
the averages are independent for different k and t, we can write 



N 



t=i 



where the averages are given as 

{dJk,t) 



exp fij + 



1 A, 



(42) 

(43) 
(44) 
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Unless we choose fij = — o"}/2, the expectation value for the jump term does not vanish. In 



the diffusion case we expanded the logarithm in Equation (11). An analogous step is not 
possible for the jump-diffusion, because the jump term dJm,t does not scale with At. This 
complicates a further analytical discussion of the jump-diffusion case. 

B. Monte Carlo results 



In the Monte Carlo simulation we consider the stochastic process in Equation (40) for 
discrete time increments. As in the diffusion case we simulate K = 500 different realizations 
of Ek^t and dJk,t for a single realization of the market terms r/m,* and dJm,t with t = 1, . . . , N. 
For each realization of the idiosynchratic terms we calculate the market return X^, the 
number of defaults A^ol-^m), the default probability PD(-^m) and the expected recovery rate 
The full Monte Carlo simulation consists of 10^ realizations of the market terms. 

For the diffusion terms and the contract details we use the same parameters as in Sec- 
tion [2j The additional parameters for the jump terms are A = 0.005, /ij = 0.4 and aj = 0.3. 

Figure |5] shows the dependence of recovery rates {R{X^)) on default probabilities PniXra) 
in one plot, and the dependence of portfolio losses (L(Xjn)) on default probabilities PdI^hi) 
in a second plot. In both cases we observe a very good agreement between the Monte Carlo 
simulations of the jump-diffusion and the analytical results for the diffusion, see Equations 



(26) and (27), respectively. The average recovery rates deviate only for very low default 
probabilities from the diffusion result. For the portfolio loss this deviation is suppressed 
and no longer visible. The parameter B is here the same as in the diffusion case, this 
means it is not influenced by the additional jump terms of the process. The most striking 
difference to the plots in Figure |2] is the abundance of simulation results with very high 
default probabilities. This is due to the market-wide jumps which render the tails of the 
market return distribution fatter. The idiosyncratic jumps lead to slightly larger deviations 
of individual results from the mean values. 

The analytical results for the diffusion case do not only describe the Pq dependence 
of the jump-diffusion results, but also their dependence on the market return X^. This 



is shown in Figure |6 



Both default probabilities and portfolio losses are well described by 



Equations (17) and (20), respectively. Given the functional relation between portfolio loss 



L and market return X^ we can again use Equation [37] to transform the pdf of market 
returns into the pdf of portfolio losses. Here we use the simulation data on X^ to determine 
their distribution numerically. In Figure [7] we compare this to the Monte Carlo result for 
the loss distribution. As in the diffusion case, we observe an excellent agreement, even for 
extremely large losses. The market-wide jumps lead to a very slow decay in the tail of the 
loss distribution. The results for Expected Loss, Value at Risk and Expected Tail Loss are 
listed in Table [Bj The analytical values describe the simulation results very well. 

Our results show that the Merton model with jump-diffusion can be described by the 
diffusion results if we exchange the distribution of market returns. The dependencies on 
default probabilities and market returns are nearly the same as in the diffusion case. This 
is somewhat reminiscent of a copula approach to credit risk. 
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Figure 5: Dependence on default probabilities Pjj for the jump diffusion process. The left plot shows the 
dependence of recovery rates {R{Xja)) on default probabilities Pu{Xjji). The right plot shows the dependence 
of portfolio losses {L{X^)) on default probabilities Py){X^). The Monte Carlo results for the jump-diffusion 
process (grey symbols) are compared to the analytical results for the diffusion (solid line). The parameter 
B = ^(1 — c)a'^T is determined from the simulation parameters. 




Figure 6: Dependence on market returns for the jump-diffusion process. The left plot shows the 
dependence of default probabilities Pd(Xi„) on market returns X^. The right plot shows the dependence 
of portfolio losses {L{Xn^)) on market returns X,„. The Monte Carlo results for the jump-diffusion process 
(grey symbols) are compared to the analytical results for the diffusion (solid line). The parameter B = 
-\/(l — cja'^T is determined from the simulation parameters. 





simulation 


analytical result 


EL 


1.03 • 10-=^ 


9.07 - 10-^ 


VaRo.99 


1.50-10-2 


1.47-10-2 


ETLo.99 


3.70 - 10-2 


3.75 - 10-2 



Table 2: Expected Loss, Value at Risk and Expected Tail Loss for the jump-diffusion process. 
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Figure 7: Probability density function pl{L) for the portfolio losses for the jump-diffusion process. The 
black solid line corresponds to the simulation results. The grey dashed line shows the analytical results, 
where the pdf of market returns is transformed. 



4. Correlated GARCH process 

Finally, we examine the Merton model with a GARCH(1,1) model for the underlying 



asset value process. The GARCH model was first introduced by BoUerslev (1986). It is able 
to reproduce many of the stylized facts found in empirical financial time series. In particular, 
it exhibits volatility clustering and fat-tailed return distributions. The GARCH(1,1) model 
is a discrete time process with an autoregressive volatility. The return at time t reads 

rk,t = (yk,t iycTlt + Vl - cek,t) 

^k,t = «o + tti rl^t-i + Pi ^k,t~i ■ (45) 

where r]t and Ek^t are independent normal distributed random variables. The parameters ao, 
ai and /3i are chosen in order to mimic the behavior of a typical empirical time series of daily 
returns. The initial values for the volatilities have been set homogeneously as (7^,0 = o""\/At, 
where a = 0.15 is the same value used in the diffusion and the jump-diffusion case. However, 
the volatility does not remain the same for all k as t evolves, since ak,t also depends on the 
idiosyncratic random part in Vk^t-i- Thus, the homogeneity of the portfolio is lost to some 
degree. Instead the process covers a wide range of volatilities both within one realization of 
the market terms, i.e., within one portfolio, and also over different market realizations. The 
GARCH process is therefore the most general case we examine here. It can provide a good 
indication of how broadly our analytical results for the diffusion case can be applied. 
For the GARCH process the market value of company k at maturity reads 

N 

VkiT) = Vol[il+fxAt + rk,t) . (46) 
t=i 



The deterministic drift term can also be included in Equation (45) instead. Then the 
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parameters ao, ai and (3i have to be adjusted accordingly. The drift constant is again set 
to fi = 0.05. 

As in the diffusion and jump-diffusion case we simulate 10^ portfolios of size K = 500. 
For each portfolio we calculate the market return X^, the number of defaults A^D(-^m), the 
default probabilities Pb{X^) and the expected recovery rate (-R(Xm)). 

In Figure |8] we present the Monte Carlo results for the Pd dependence of the GARCH 
process. One plot shows the dependence of recovery rates (P(Xin)) on default probabilities 
PD(Xm), a second plot shows the dependence of portfolio losses (L(Xm)) on default prob- 
abilities PD(Xm). First, we observe a much broader range of the individual results than in 
the diffusion and jump-diffusion case. This is to be expected because the fluctuating volatil- 
ities lead to pronounced fat tails in the return distribution of individual companies. While 
the returns are correlated in the GARCH(1,1) model, volatilities are rather uncorrelated, 
see Schafer and Guhr (2010). Thus, for the market return the fluctuations of individual 
volatilities average out to some degree. This is why we scarcely observe very high default 
probabilities in the simulations. The situation is different for empirical stock return time 
series, where volatilities are also strongly correlated. 

It is a rather remarkable result that the average behavior of the GARCH simulation is 



so well described by the analytical results for the diffusion, see Equations (26) and (27), 
respectively. The parameter B has been fitted to the simulation data. 

In Figure [9] we demonstrate the dependence of default probabilities and portfolio losses 
on the market return. Here the behavior of the GARCH model is completely different from 
the diffusion case. For the analytical results of the diffusion case, described by Equations 



(17) and (20), we use the parameter value B which has been fitted to the Pd dependence in 
Figure [8j Given the functional relation between portfolio loss L and default probability Pd 
we can transform the pdf of default probabilities into the pdf of portfolio losses, 

Pl(I^) = ^^7(^J>Po™. (47) 
A comparison to the Monte Carlo result for the portfolio loss distribution is presented in 



Figure [TOj We observe a very good agreement, even for extremely large losses. Finally, we 
present the results for Expected Loss, Value at Risk and Expected Tail Loss in Table |4} We 
compare the Monte Carlo result for Value at Risk with 

VaR, = L(Pp^i(a)), (48) 

where Fp^{a) is the a-quantile of default probabilities. The Expected Tail Loss can be 
expressed as 

ETL, = - / L(PD)pp^(PD)dPD . (49) 
a J 

The analytical values describe the simulation results very well. 

Our results demonstrate that the functional dependence of recovery rates on default 
probabilities does not depend on the underlying process. The analytical result for the 
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Figure 8: Dependence on default probabilities Pd for the GARCH process. The left plot shows the 
dependence of recovery rates {R{X^)) on default probabilities PuiX^). The right plot shows the dependence 
of portfolio losses {L{X^)) on default probabilities Py){X^j^). The Monte Carlo results for the GARCH process 
(grey symbols) are compared to the analytical results for the diffusion (solid line). The parameter B is fitted 
to the data. 



diffusion process describes also other underlying processes in the Merton model. Hence we 



call Equation (26) the structural recovery rate. Using this structural recovery rate we can 



describe the loss distribution if we have knowledge of the distribution of default probabilities. 





simulation 


analytical result 


EL 


2.44 ■ 10-3 


2.26 • 10-3 


VaRo.99 


3.28 • 10-2 


3.16 ■ 10-2 


ETLo.99 


5.42 ■ 10-2 


5.16 ■ 10-2 



Table 3: Expected Loss, Value at Risk and Expected Tail Loss for the GARCH(1,1) process. 



5. Applications of the structural recovery rate 



In a recent study Altman et al. (2005) find a strong negative correlation between default 
probabilities and recovery rates in empirical data. The structural recovery rate provides 
such a negative relation between default and recovery rates. It correctly describes the case 
of zero-coupon bonds for various underlying processes, but may even be apphcable to a 



more general debt structure. An indication for this is given by Chen and Panjer (2003) who 



show that the Merton model with jump-diffusion can be calibrated to provide the same 
yield spreads as reduced-form models. 

Reduced-form models are also able to reproduce negative correlations between default 
and recovery rates, if the default model and the recovery model depend on a single common 
covariate, see Chava et al. (2008). Compared to the structural recovery rate, however, the 



reduced-form recovery introduces more parameters and lacks a deeper motivation. 
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Figure 9: Dependence on market returns for the GARCH process. The left plot shows the dependence of 
default probabilities PuiXj^^) on market returns The right plot shows the dependence of portfolio losses 
on market returns X^-i- The Monte Carlo results for the jump-diffusion process (grey symbols) are 
compared to the analytical results for the diffusion (solid line). The parameter B is fitted to the default 
rate dependence in Fig. |8] 




0.00 0.05 0.10 0.15 0.20 

L 

Figure 10: Probability density function pl{L) for the portfolio losses for the GARCH process. The black 
solid line corresponds to the simulation results. The grey dashed line shows the analytical results, where 
the pdf of default probabilities is transformed. 
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Figure 11: Empirical data from Moody's annual default study for the years 1987 to 2008. The left plot 
shows the dependence of firm-wide ultimate recovery rates on the annual default rates. The right plot shows 
the dependence of the corresponding average losses on the annual default rates. The solid line shows the 
structural recovery rate with the parameter B = 2.28 fitted to the data. 



In first passage models as described in Giesecke (2004), default may occur before maturity 
if the market value of the company falls below a default barrier. If defaults before maturity 
dominate in this model, recovery rates are independent from default rates. The average 
recovery rate is then constant and does not resemble the negative relation mentioned above. 
However, if defaults occur mostly at maturity, the structural recovery rate is recovered. 

Finally, we give a first empirical indication that the structural recovery rate is indeed 
applicable in a realistic setting. Figure 11 shows data from Moody's annual default study. 
Although there are only very few data points, corresponding to the years 1987 to 2008, 
we can observe a negative relation between the firm-wide ultimate recovery rates and the 
annual default rates. We fit the structural recovery rate to these data points and find a 
reasonable agreement. 



6. Conclusions 

The interdependence of default and recovery rates has a crucial infiuence on large credit 
losses. Yet default probabilities and recovery rates are often modelled independently in cur- 
rent credit risk models. In this paper we revisited the Merton model for different underlying 
processes with correlations. While the original Merton model was conceived for the diffusion 
process, we also consider a jump-diffusion and a GARCH process. For the correlated diffu- 
sion we derived a functional dependence between default and recovery rates. This functional 
dependence is determined by a single parameter. In Monte Carlo simulations we showed 
that it describes not only the diffusion case, but also the jump-diffusion and GARCH pro- 
cess. Due to its independence on the underlying process, we call this functional relation the 
structural recovery rate. It is straightforward to use this structural recovery rate in addition 
to a first passage model or any other model for default probabilities. We believe that this 
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has great potential to improve current credit risk models. In further studies we shall address 
two important questions: How well does the structural recovery rate describe empirical data 
on defaults and recoveries? And how does it compare with reduced-form models? 
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